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^ ; ABSTRACT 

Context. The use of model atmospheres for deriving stellar fundamental parameters, such as Teff , log g and [Fe/H] , will increase as we 
' find and explore extreme stellar populations where empirical calibrations are not yet available. Moreover calibrations for upcoming 

I large satellite missions of new spectrophotometric indices, similar to the uvby-H/} system, will be needed. 

^ , Aims. We aim to test the power of theoretical calibrations based on a new generation of MARCS models by comparisons with 

observational photomteric data. 

^ Methods. We calculate synthetic uvby—\ip colour indices from synthetic spectra. A sample of 388 field stars as well as stars in 

I globular clusters is used for a direct comparison of the synthetic indices versus empirical data and for scrutinizing the possibilities of 

theoretical calibrations for temperature, metallicity and gravity. 
I ■ Results. We show that the temperature sensitivity of the synthetic (b - y) colour is very close to its empirical counterpart, whereas 

O ' the temperature scale based upon Hj8 shows a slight offset. The theoretical metallicity sensitivity of the mi index (and for G-type stars 

^ ' its combination with ci) is somewhat larger than the empirical one, based upon spectroscopic determinations. The gravity sensitivity 

C/3 , of the synthetic C[ index shows a satisfactory behaviour when compared to obervations of F stars. For stars cooler than the sun a 

deviation is significant in the Ci-{b - y) diagram. The theoretical calibrations of (b - y), (v - y) and ci seem to work well for Pop II 
' stars and lead to effective temperatures for globular cluster stars supporting recent claims by Korn et al. (2007) that atomic diffusion 

^ ' occurs in stars near the turnoff point of NGC 6397. 

, Conclusions. Synthetic colours of stellar atmospheres can indeed be used, in many cases, to derive reliable fundamental stellar 

. parameters. The deviations seen when compared to observational data could be due to incomplete linelists but are possibly also due 

' to effects of assuming plane-parallell or spherical geometry and LTE. 

CO ' Key words, stars: atmospheres - stars: synthetic spectra ~ stars: fundamental parameters - techniques: photometric 

O ; 
(N ' 

1. Introduction poor stars. In the calibration by Alonso et al. (1996) the effec- 
QP tive temperature calibration was based on stellar temperatures 
O The uvby photometric intermediate-band system of Stromgren fj-om the Infrared Flux Method. The calibrations by Nordstrom 
(1963) and the H/3 narrow-band system (Crawford 1958, gt al. (2004) and Holmberg et al. (2007) additionally utilised 
K~Ji ■ ^^^^^ ^^^"^ combined early and soon proved to be a most Hipparcos parallaxes for the surface gravity calibration. 
pX^ powerful means of determining fundamental parameters of Another development towards calibration of the systems also 
H stars. Calibrations of the systems were also developed early gj^rted early: the direct calculation of photometric indices by 
. F^. by Stromgren, Crawford and collaborators. Generally, semi- means of model atmospheres. Such theoretical calibrations were 
empirical methods were used with the calibrations made by attempted for early-type stars with relatively line-free specti-a. 
means of sets of stars with parameters determined in more fun- p^^. i^jg jypg ^^^.^^ ^ statistical correction for the effects of spec- 
damental ways. The first metallicity calibration of the mj index ^.^^ jings was made by Baschek (1960) in his calibration of the 
for solar-type stars by Stromgren (1964) was thus based on spec- stromgren m index (a predecessor to mi). A first systematic and 
ti-oscopical [Fe/H] values of Wallerstein (1962) and the first lu- detailed calculation of uvby-U/3 indices for a grid of F and G 
minosity calibration of the ci index by Stromgren (see Crawford ^jwarf model atmospheres was published by Bell (1970), using 
1966) used cluster stars in the main-sequence band. This general gg^led solar model atmospheres. Bell & Parsons (1974) calcu- 
semiempirical approach continued with calibrations like those i^tej ^^^^^^^.^ for flux-constant model atmospheres of F and 
of the mi index by Nissen (1970), Gustafsson & Nissen (1972), q supergiants, while Gustafsson & Bell (1979) produced theo- 
Nissen & Gustafsson (1978) and Nissen (1981) where very- ^etical colours in a number of systems, including the uvby sys- 
naii-ow-band spectrophotometi-y of groups of weak lines was ^^^^ fo^ a grid of giant-star model atmospheres. Relya & Kurucz 
cahbrated with synthetic spectra. These stars were then used as (1973) calculated uvby and UBV colours from early ATLAS 
calibration stai-s for Mvfey-H^ photometry. Other calibrations fol- models, and discussed their shortcomings for late-type stars, 
low this semiempirical approach e.g. those of Crawford (1975), ^^^,l,y colours for new sets of Kurucz models were published by 
Ardeberg & Lindgren (1981) and Olsen (1988) for G and K- Lester, Gray & Kurucz (1986). Castelli & Kurucz (2006) pub- 
type stars, as well as by Schuster & Nissen (1989) for metal- ijghed H/3 indices. Sometimes, semiempirically corrected fluxes 

from model atmospheres have also been used for calibrations of 
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The need for reliable calibrations of uvby-UJi photometry 
has increased in the last decade, not the least for estimating pa- 
rameters of new and more "exotic" stars, such as very metal-poor 
and super-metal-rich stars which are not found at great abun- 
dance in the solar neigbourhood so that relatively complete sets 
of calibration stars cannot easily be established. Furthermore, 
the preparation for the Gaia satellite includes a careful analysis 
of the power of model atmospheres to provide a detailed astro- 
physical calibration of the photometric system of the satellite. 
This analysis needs support by a detailed test of the problems 
and possibilities to make a detailed theoretical calibration of, 
e.g., the uvby-li/5 photometry. 

Subsequently, we shall present the theoretical models and 
colours (Sect. |2] and O. Stellar samples for empirical compar- 
isons are discussed in Section |4] Next, the discussion will be 
focused on the determination of effective temperature, metallic- 
ity and surface gravity of the stars, by discussing the calibration 
of (b - y) and H/3 indices (Sect.|5] effective temperature), of mi 
(Sect. |6] metallicity) and ci (Sect. |7] gravity) indices, devoting 
more limited interest to "secondary" effects such as the metallic- 
ity sensitivity of (b - y) and ci, or the gravity sensitivity of nii. 
In each section the results will be compared with empirical and 
semi-empirical data and calibrations. Finally, in the last Section 
some comments will be made on the success and the problems of 
the theoretical calibrations, conclusions will be drawn and rec- 
ommendations given. 



2. Model atmospheres and calculated spectra 

The theoretical tools used in modelling the stellar colours are 
model stellar atmospheres and their calculated fluxes. These are 
based on extensive atomic and molecular data. Here, we shall 
briefly present the models and data used and refer to more com- 
plete descriptions. 



2.1. Model atmospheres 

The stellar atmosphere code MARCS (Gustafsson et al. 2008, 
\http://marcs. astro. uu.se) was used to construct a grid of 168 theo- 
retical ID, flux constant, radiative + mixing-length convection, 
LTE models with fundamental atmospheric parameters as fol- 
lows: Teff = 4500, 5000, 5500, 6000, 6500 & 7000 K, \ogg = 
2.0, 3.0, 4.0, 4.5 and [Me/H] = 0.5, 0.0, -0.50, -1.0, -2.0, - 
3.0, -5.0, [Me/H] denoting the logarithmic over-all metallicity 
with respect to the sun. Plane parallell (ppl) stratification was 
assumed for log g - 4.5, 4.0 & 3.0, whereas spherical (sph) sym- 
metry was assumed for logg =2.0. For the spherically symmet- 
ric models a mass of 1 M© was adopted. The local mixing-length 
recipe was used to describe convective fluxes, for more details 
see Gustafsson et al. (2008). 

Elemental abundances were adopted from Grevesse and 
Sauval (1998) except for the CNO abundances which were 
adopted following Asplund et al. (2005). 



2.2. Synthetic spectrum calculations 

In order to calculate synthetic spectra of sufficiently high resolu- 
tion the Uppsala BSYN code was used with the MARCS models 
as input. The spectra were calculated within the wavlength limits 
of the Stromgren uvby filters with wavelength steps of 0.02 A. 
The mictroturbulence parameter, was set to l.Okms"' and 
1.7kms"\ when calculating spectra based on ppl models and 



sph models, respectively. The effects of changing ^, are studied 
in Section|6] 



2.3. Line lists 

We collected atomic line data from the Vienna Atomic Line 
Data Base, VALD (version I, Kupka et al. 1999). For hydro- 
gen line data a version of the code HLINOP was used, and has 
been described by Barklem & Piskunov (2003). This code has 
been developed based on the original HLINOP by Peterson & 
Kurucz (see http:// kurucz.harvard.edu/). The hydrogen line pro- 
files are calculated including Stark broadening, self-broadening, 
fine structure, radiative broadening, and Doppler broadening 
(both thermal and turbulent). The Stark broadening is calculated 
using the theory of Griem (1960 and subsequent papers) with 
corrections based on Vidal et al. (1973). Self-broadening is in- 
cluded following Barklem et al (2000) for Ha, li/3 and Hy, while 
for other lines the resonance broadening theory of Ali & Griem 
(1966) is used. For molecules, data for C2 including '^C'"'C lines 
were gathered from Querci, Querci & Kunde (1971) and Querci 
(1998, private communication with B. Plez), except for the Fox- 
Herzberg band in the UV for which data was taken from Kurucz 
(1995). Molecular data for CH with '^CH comes from Plez et al. 
(2008) and Plez (2007, priv. communication). Data for the CN 
molecule with '^C'^N, '^C"*N, and '^C'^N are also from Plez 
(priv. communication). CO data with '^C'^O was gathered from 
Kurucz (1995) as well as NH with ^^NH, OH with "^OH, and 
MgH with ^^MgH and ^^MgH. TiO was not taken into account in 
the calculations since absorption by this molecule is expected to 
have small effects on the spectra within the effective-temperature 
range considered. 



3. Colour index calculations 

3.1. Filter profiles 

To determine the theoretical mi, ci and (b-y) indices, transmis- 
sion profiles of the Stromgren uvby filters (Crawford & Barnes 
(1970) see Fig. [U were multiplied with the calculated model 
stellar surface flux within the wavelength range of the filters: 



mag = —2.5 log 



( J F,T,dA 
J TAdA ) 



+ const. 



where and T,i are the flux and the relevant transmission 
profile, respectively. The theoretical magnitudes were converted 
into colour indices via the definitions: ci = u - 2v + b, mi = 
V — 2b + y. 

The H/3 index is defined (Crawford 1958) as the ratio of the 
flux measured through a narrow and a wide profile, respectively, 
both centered around the H/3 line: 



H/3 = -2.5 



log log 



/ TN,AdA 



J TwjdA ^ 



+ const. 



where Tw,a and Tn,a are the wide and the narrow filter profiles, 
respectively. We have found notable differences between syn- 
thetic Hj6 indices calculated with the Kitt Peak filter set denoted 
(9,10) and those calculated with the (212,214) set, both de- 
scribed by Crawford & Mander (1966, CM66). Here, the Kitt 
Peak (212,214) filter system was chosen (Fig. 1 CM66) with a 
transmission function adopted from CastelU & Kurucz (2006), 
see Figure [U 
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Fig. 1. The uvby-li/3 transmission functions of the standard sys- 
tems plotted as a function of wavelength. As a comparison, 
the flux (per Angstrom unit) of a model with Teff = 6000 K, 
log g - 4.0 and [Me/H] - 0.0 is plotted on an arbitrary flux scale. 



3.2. Transformation to the observational system 

The frequent use of filter-profiles diff'erent from those that 
originally have defined the photomteric system needs some 
extra considerations. For the H/3 index this mainly comes in at 
the zero-point determination. Since the observed indices are 
defined on a system with zero-points set by particular standard 
stars, we must apply corresponding zero-point shifts to our 
theoretical system. The much observed, and frequently used, 
star Vega was chosen to estimate these zero-points. Thus, a 
model atmosphere for Vega was calculated with parameters Teff 
= 9550 K, logg = 3.95 and [Me/H] = -0.5 (mean values of 
selected measurements presented in SIMBAD), and a synthetic 
spectrum was computed. Theoretical indices for the Vega model 
were calculated and compared to observed values: ci - 1.088, 
mi =0.157 and (b - y) = 0.003 (Hauck and MermilHod 1998). 
The resulting differences between the theoretical and observed 
colours were then added as constants to all calculated Stromgren 
colours in our model grid. 

For the H/S index, Crawford & Mander (1966) presented several 
filter systems, of which we choose the (212,214) filters as 
mentioned above. Indices calculated by using this filter set, 
which are referred to as H"P below, should be transformed 
via a set of equations (CM66), in order to agree with previous 
Crawford-Mander observations on their standard system. The 
transformation for the (212,214) filters is described by the 
fonowing two equations (CM66, Table III): 

np = 0.374 + 1 .305 H'/3, B stars 

H/3 = 0.248 + 1 .368 H'/3, A, F stai's, 

derived from a set of 45 and 35 bright stars, where H'/3 
and Hy6 are the observed and tranformed indices, respectively. 
The later equation for A and F stars was applied on all H/3 
indices in our theoretical grid. 

The AO star Vega clearly poses a number of problems for 
determining the zero-point of the uvby and H/3 indices. It is 
known to be rapidly rotating, but with its axis close to the line of 
sight (Gulliver, HiU & Adelman 1994, Hill, Gulliver & Adelman 
2004). Vega has also been regarded to show mild A Bootis star 



characteristics, such as certain non-solar abundance ratios as 
well as dust emission in the IR (see Gigas 1988, Hill 1995, 
Ilijic et al. 1998, Adelman & Guflliver 1990, Heiter, Weiss 
& Paunzen 2002). However, these departures from standard 
AO stars, as well as standard model atmospheres, are thought 
to only lead to minor modifications of its uvby-Hp indices 
(Paunzen et al. 2002). A more practical problem is that neither 
of the two CM66 transformation equations from H'/J to H/3 
will give a fully satisfactory fit due to the fact that the index 
for Vega should be transformed intermediately between the B 
star and the A, F star sets. Several tests were performed which 
all pointed in the direction that Vega should be transformed to 
the standard system via an equation somewhere in between the 
B and the A, F transformations but with a heavier weight for 
the latter By using all Usted AO stars in Table II of CM66, a 
special transformation between H'/3 and H/3 for AO dwarf stars 
was established, as can be seen in Figure |2] The vertical line 
in the plot represents the observed H/3 value for Vega (Hauck 
& Mermilliod, 1998). However, the observed H'/3 value of 
Vega is not known. In order to determine the zero point c in 
the transformation from calculated {H"/3) to "observed" H'/3 
values, H"/3 = H'/3 + c, we have therefore adopted an H'yS value 
for Vega derived from the AO stars line in Figure |2] read ofi" at 
the observed Hy6 for the star. After correction of the model H"y6 
value to H'y6 we have then calculated the model HyS values by 
using the empirical A,F-transformation relation as mentioned 
earlier 



B stor equotion 
A,F stor equotion 
AO star linear fi: 




2,85 2.88 



Fig. 2. H/3 and H'/3 values for AO stars tabulated in CM66 are 
plotted individually and represented by a linear regression (thick 
line), as well as the transformation relations for B and A+F stars, 
respectively, following CM66 (thin lines). 



The transmission functions of the filters used by Olsen 
(1983, 1984) and Schuster & Nissen (1988) depart from those on 
which the system was originally based (Crawford 1966). Filter 
profiles representative for these newer studies are given by Helt 
et al. (1987) and Bessell (2005). In particular, the more recent v 
band is narrower, and its effective wavelength shifted by about 
25 A towards the red. To test the eflects of this, we have calcu- 
lated colours with the Helt et al. (1987) profiles as an alterna- 
tive to the Crawford set, and then applied the transformations 
as given by Schuster & Nissen (1988) back to the standard sys- 
tem to mimic the procedure of the observers. From this, we find 
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changes of the calculated values, amounting to typically 0.02 
magnitudes or less in the mi and ci index values and only 0.002 
in (b - y). In particular, the changes in the differential values 
measuring the sensitivity of mi to metallicity and of ci to grav- 
ity are small, in 6mi / 6[MefH] typically 0.003 and in 6ci/6logg 
typically 0.01. Such effects do not change our conclusions in the 
present paper. 

3.3. Model colours 

The model colours, with zero-point added using Vega observa- 
tions, are supplemented the present paper electronically. 



4. Comparison star samples 

In order to test the reliability of our calculated colours, we se- 
lected a sample of standard stars. These were taken from vari- 
ous sources: one subset with well determined spectroscopic pa- 
rameters was selected from The Bright Star Catalogue (Hoffleit 
& Warren 1995), another from the uvby standard stars listed 
by Crawford & Barnes (1970), a third from the list of metal- 
poor stars in Schuster & Nissen (1988), a fourth among stars 
that have been observed by the Hubble-STIS spectrograph and 
a fifth from the stars listed in the study of Pop II stars by 
Jonsell et al. (2005). Altogether 388 stars were thus selected. 
The fundamental parameters were taken from the sources listed, 
or from other sources given in the SIMBAD catalogue and 
judged to have high quality. Complementary photometry was 
also obtained from SIMBAD. Parameter determinations based 
on uvby-H/3 photometry were avoided as far as possible, since 
we aimed at testing calibrations based on this photometric sys- 
tem relative to parameters based on more fundamental methods. 
In practice, this usually means effective temperatures based on 
the infrared-flux method and gravities and metallicities based on 
high-resolution spectroscopy. The effective temperatures gath- 
ered from the Jonsell et al (2005) sample are calculated with the 
Alonso et al. (1996) calibration. These stars however, constitute 
less than 5 % of our total sample. In the case of multiple sources, 
i.e., stellar values listed in more than one of our selected cata- 
logues, a mean value was used. The uvby indices of the standard 
stars were dereddened by means of the algorithm and computer 
code of Hakkila (1997). For most of the stars, in particularly the 
dwarfs, ^(B-V) was smaller than ~0.01. 

Altogether, the standard stars span a volume in the param- 
eter space ranging from 3900 K to 7850 K in Teft, 0.20 to 4.80 
in \ogg (cgs units) and -3.0 dex to 0.45 dex in [Fe/H]. Data for 
the full standard sample is accessible electronically from A&A 
as on-line material supplementing the present paper As seen in 
Figure |3] our standard sample satisfactorily covers the funda- 
mental parameter space, although one would need some more 
cool (reft < 5000) and hot (T^ff > 6300) dwarf stars (log^ > 
4.0). For this purpose two complementary homogeneous sam- 
ples were tested; the stellar sample published by Casagrande 
et al (2006, hereafter C06) and Valenti & Fischer (2005, here- 
after VF05). uvby colour indices for both samples were taken 
from Hauck & Mermilliod (1998). These complimentary stel- 
lar samples cover the parameter spaces: 4700 < Tf^g < 6600, 
3.1 < log ^ < 5.1 and -1.9 < [Fe/H] < 0.5 for the VF05 sample 
and 4406 < T^g < 6556, logg ~ 4.5 and -1.87 < [Fe/H] < 0.34 
for the C06 sample. Despite that the lowest logg values in the 
VF05 sample are characteristic of giants/sub-giants, a majority 
(93 %) of the VF05 stars are dwarfs with log g > 4.0. In the C06 
sample all stars are assumed to be dwarfs with logg ~ 4.5. A 



majority (85 %) of the stars in the VF05 sample are also metal- 
rich, [Fe/H] > -0.2, whereas the C06 sample is evenly dis- 
tributed in metallicity. Likewise, the effective temperatures of the 
C06 sample (determined by the authors's new calibration of the 
infrared-flux method) are evenly distributed but the VF05 sam- 
ple is biased towards higher temperatures, T^ff > 5500 (74 %). 
Among the stars in the standard sample, 18 and 72 are found 
in the C06 and VF05 samples, respectively. For the stars in the 
VF05 sample no significant trends of discrepacies for the val- 
ues given, Tf.ff, logg and [Fe/H], were found. In the C06 sam- 
ple we find an overall difference in given effective temperatures, 
growing for higher effective temperatures at low metallicities 

(T'eff standard " T'eff Casagrande ~ -175 K at Jeff ~ 6000 K and -2.0 < 

[Fe/H] < -1.0). 




7000 6000 5000 4000 



Fig. 3. Fundamental parameters, i.e. Teff, [Fe/H] and logg, for 
the standard sample. 



5. The effective-temperature calibration 

Precise determinations of effective temperatures of stars are 
critical for a variety of reasons, not only for direct applica- 
tions such as comparisons of isochrone calculations to observed 
colour-magnitude diagrams, but also in more indirect applica- 
tions e.g. determinations of chemical abundances. Photometric 
indices play a critical role when determining the effective tem- 
peratures of stars. For stars located on the main sequence, model 
atmosphere calibrations of these indices may be particularly im- 
portant since determinations of diameters are generally few and 
poor 

The Stromgren (b-y) index provides a sensitive temperature 
measure for F, G and K stars. The H/S index is another frequently 
used criterion for deriving temperatures. We study both indices 
here. 

In Figure |4] (left) we explore the temperature sensitivity of 
the theoretical (b-y) index for different metallicities. As a com- 
parison we also plot (b - y) for a set of Kurucz models (Lester, 
Gray & Kurucz, 1986). We note that the indices, calculated with 
the different model sets, more or less show the same behaviour 
although the models of Lester et al. are systematically bluer at 
given Tfiff. In comparison with observed values for stars our syn- 
thetic indices show a satisfactory temperature sensitivity over the 
full temperature range (4500-7000 K), see Figure|5] 

In Figure |4] (right) the H/3 index is plotted versus effective 
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temperature, together with the theoretical data of Castelli & 
Kurucz (2006). For the warmer part of the temperature range 
we find no differences in sensitivity, whereas for the cooler tem- 
peratures (5500 K) we note that MARCS indices have a steeper 
gradient, thus lower H/3 values than those of Castelli & Kurucz. 
Our steeper gradient tends to agree with observed stellar colours. 
Figure |6] which is probably due to improvements of the broad- 
ening theory of H/3. 

5.1. T,nib-y) 

Several attempts (Alonso et al. 1996, Nordstrom et al. 2004, 
Ramirez & Melendez 2005, Holmberg et al. 2007) have been 
made to derive empirical calibrations for uvby-Hfi colours, 
where the fundamental parameters are expressed as explicit for- 
mulae in the photometric indices. We will here derive such cali- 
bration equations by using theoretical colours. We make use of 
the formal expressions of Alonso et al. (1996, 1999, hereafter 
A96 resp. A99) based on (b-y) (their eq. 9 and 14 & 15, respec- 
tively) but derive new coefficients via an iterative least squares 
method and by using the model colours. 

The result for dwarfs (logg > 4.0) are given numerically in 
Appendix lA. 1 l and can be seen in Figure Q where we also merge 
our standard sample with the C06 sample and plot temperatures 
for the individual stars derived from our theoretical calibration 
equation. The standard deviations for our calculated tempera- 
tures with respect to literature values can be found in Table [T] 
A linear regression for the calculated temperatures of the two 
merged samples is also shown in Figure|2]and yields a dispersion 
cr = 1 17 K. We note that our theoretical calibration gives higher 
temperatures, as compared with literature values, by ~100K for 
the cooler stars (Teff < 5000 K) and tends to lead to lower temper- 
atures (~100-150K) than the values found in the literature for 
hotter stars (Teff > 6500 K). Our calibration suggests lower tem- 
peratures than the values listed in VF05 by 1 00-200 K for stars 
within the range of 5200-6300 K. Below that our theoretical cali- 
bration gives systematically ~100K higher temperatures for the 
VF05 stars, in accordance with the tendencies found also for the 
other comparison samples. The standard deviation of the calcu- 
lated temperatures for the VF05 sample can be seen in Table[T] 

It is also interesting to compare our results to previous empir- 
ical calibrations. The standard deviations of the calculated tem- 
peratures using four different empirical calibrations applied on 
our three comparison samples, are listed in Table [T] In order to 
illustrate the empirical trends the resultant linear regressions for 
the merged standard and C06 samples are shown in Figure |7] 
(note the shifted scale by 1000 K). 

The calibration results for giants (1.5 < \ogg < 3.5) is pre- 
sented in Appendix IA.2I and shown in Figure [8] The standard 
deviations with respect to effective temperatures of the standard 
stars are shown in Table [T] together with the standard deviations 
for the same sample when using the A99 empirical calibration. 

5.2. T,fim 

A theoretical temperature calibration for dwarfs was also derived 
for the H/3 index, with a formal expression adopted from A96, 
Eq. 10, see Appendix lA.il The result can be seen in Figure |9] 
Standard deviations with respect to the literature values of the 
standard sample are listed in Table [1] values somewhat higher 
than the calibration based on (b-y). Here we are, however, using 
different stars in the T^^ffib - y) and Teff (HjS) equations (fit limits 
defined in A96). To make a fair comparison between the two 



Table 1. Standard deviations for Tf.fi and [Fe/H] derived using 
models as compared with literature values of the different com- 
parison samples. 
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calibrations we have restricted our sets of stars to be identical 
for both calibrations. Such a comparison is made in Sec. 15.31 
A linear regression for the calculated effective temperatures is 
also shown in Figure |9] The dispersion around this fit is cr = 
143 K. We see that our theoretical calibration suggests somewhat 
higher effective temperatures, of roughly lOOK, in the cooler 
part of the temperature range (Teff ~ 5000 K), and that it possibly 
implies lower effective temperatures of the order of ~ lOOK for 
the warmest stars (Tgff > 6500 K). Some of these tendencies may 
possibly be traced also in the (b - y) - T^fi calibration (c.f. Fig. 
|7]i. As an alternative, the effective temperatures for the standard 
sample are calculated with the empirical calibration equation by 
A06. The standard deviations with respect to literature values 
are presented in Table[T] and a linear regression to the calculated 
effective temperatures is also plotted in Figure |9] 

5.3. T,ff(b-y)vsT,ff(H/3) 

Now we compare the derived temperatures based on (b - y) and 
H/3, respectively. 119 stars matching the restrictions in the pa- 
rameter space for both calibrations were selected out of our stan- 
dard sample. The effective temperatures based on (b-y) and Hy6 
and the result are displayed in Figure [10] The test reveals that 
there is some disagreement between the two theoretical calibra- 
tions. The H/3 based equation suggests higher temperatures for 
the cooler stars (Teff < 5000 K) and lower temperatures for the 
warmer stars (Teij < 6900 K) of some 200 K. The empirical equa- 
tion of Alonso et al. (1996) shows more or less the same spread 
for the standard sample but no significant departures from the 
1-1 line, which is also to be expected. The failure for the theo- 
retical calibration for the coolest stars is not very remarkable, 
since the metal lines for those stars are strong and dominate the 
H/3 line, and other parameters affect the index such as gravity. 
The departure in the hotter end is only dependent on a few stan- 
dard stars. 

It is worth noting that the difference between the (b - y) 
and HyS calibrations responds in opposite directions to differ- 



6 



Onehag et al.: Calibration of Stromgren Piiotomtery 



0.8 
0.6 
0.4 

0.2 

0.0 
0.8 

0.6 

0.4 

0.2 

0.0 
0.8 

0.6 

0.4 

0.2 
0.0 



[K/e/H] = 0.00 ^ : MARCS 

X X : LGK 




4500 5000 5500 6000 6500 7000 4500 5000 5500 6000 6500 7000 



Fig. 4. The (b - y) colour index (left) and H/3 index (right) versus effective temperature for MARCS and Kurucz (Lester, Gray & 
Kurucz (1986), LGK, Castelli & Kurucz, 2006, CK) model atmospheres with different [Me/H] and \ogg (indicated to the left and 
right of each curve). Corresponding theoretical MARCS {b - y) and H/3 curves for [Me/H] = 0.5, -0.5, and -3.0 can be seen in 
Figure |5] and Figure |6] respectivley. 



ent gravities. Le., the temperature sensitivity 6{index)/S(Tf.f[) 
decreases with increasing gravity for {b - y) while it increases 
for Teff < 5000/: for ¥1/3 (see Fig. gli. As is seen, these differ- 
ences are also metallicity dependent, and tend to vanish for small 
metallicities. However, for stars with known reddening and as- 
suming e.g. the metallicity to be known and relatively large, this 
may make it possible to obtain a temperature and rough gravity 
classification from (b - y) and HyS, only. 

6. The metallicity calibration 

The basic aim of nii = (v - b) - (b - y) was to measure the 
total intensity of the metal lines in the v band. As was early 
appreciated by Stromgren, for late F and G stars of Pop I these 
lines are, however, to a large extent located on the flat part of the 
Curve of Growth and are thus not very sensitive to metallicity, 
but rather to micro turbulence. Moreover, for the hotter stars, the 
US line is strongly affecting the band, and for stars later than 
G5, CN lines of the (0,1) band in the Violet System are also 
significant. Thus, the effects of the value of the mircroturbulence 
parameter, as well as of the individual CNO abundances, e.g. 
due to dredge-up of CNO processed material from the interior, 
must be taken into consideration. Certainly, mi also varies with 



effective temperature and, to a less degree, with surface gravity. 

The variation of the calculated mi with (b - y) and metal- 
licity for the model atmospheres is shown in Figure [TT| and 
compared with calculated indices by Lester et al. (1986). As is 
seen, the index offers a good discrimination in metallicity ex- 
cept for stars of Extreme Population II for which it only works 
for the cooler end of the temperature interval. A characteristic 
measure of the sensitivity of the index to overall metallicity is 
{6{mi) / 6[Me/H])^b-v), where the subscript denotes that the sen- 
sitivity is measured at a constant (b - y). This quantity, as mea- 
sured for models with -3.0 < [Me/H] < 0.5 and logg ^4.5 & 
4.0, is given in Table |2] 

We may compare the sensitivity of the mi index with the 
empirical results. However, recent calibrations of the Stromgen 
photometry, e.g. by Holmberg et al. (2007), contain complex 
non-linear expressions in which all the indices are involved - 
in principle a reasonable approach since e.g. also the ci index 
carries information on the metallicity. Since this metallicity de- 
pendence is, however, far from independent of that carried by 
the mi index the terms in the calibration expressions involving ci 
may well mask some of the dependence of mi on metallicity. So, 
as we wish to understand the way mi changes with metallicity. 
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Fig. 5. Theoretical (b-y) values versus effective temperature for models with different [Me/H] and log § plotted together with values 
for the standard stars (with [Fe/H] within +0.1 dex of the given [Me/H]). To guide the eye the theoretical log^ curves for 2.0 and 
4.5 are plotted as a dashed and a dotted line, respectively. 



we have instead turned back to the earlier empirical caUbrations 
like that of Nissen (1988), where the metallicity dependence of 
mi was still treated separately. 
Thus, Nissen (1988) finds 

[Fe/H] = -(10.5 + 50(HyS- 2.626)) -too + 0.12 (1) 

where dniQ is calculated as nistand - '"i relative to the Hyades 
standard sequence at constant HyS and corrected for interstel- 
lar extinction. From Figure [TT] or Table |2] we may estimate 
5OTi/(5[Me/H] at constant {b - y). We can easily calculate the 
corresponding empirical quantity from the coefficient in Eq. ([T]i 
(correcting for the derivative at constant HyS to constant (b - y) 
by adding 0.03 which is easily found with sufficient accuracy 
from the model colours) and then find empmcal sensitivities 
5mi/5[Fe/H] of 0.093 (7000 K) and 0.15 (5750 K), with rele- 
vant effective temperatures within parentheses. These values for 
the sensitivities are valid for [Fe/H]* 0.0. The corresponding 
theoretical sensitivities 5mi/5[Me/H] for dwarf stars are typi- 
cally 0.10 and 0.14, respectively. The agreement between the 
empirical and theoretical results is quite satisfactory. It seems 
possible that a basic reason why the theoretical calibration, in 
spite of this agreement, does not succeed very well for the more 
metal-poor stars (see below) may rather be due to the modelling 
of the fluxes of those stars than due to failures in calculating the 



change of flux with metallicity for solar-type stars. It may also be 
associated with the measured v band being different to that com- 
puted, as a single transformation equation is unlikely to correct 
solar metallicity stars and metal-poor stars by the same amount. 

In order to further explore the properties of the calculated 
m\ indices we have plotted individual stars with fairly well- 
determined fundamental parameters, chosen from our standard 
sample in the m\ - (b - y) diagram. As is seen in Figure [12] 
these stars match the calculated indices relatively well, although 
there seems to be a tendency of the sensitivity of the mi index to 
metallicity to be exaggerated by the model fluxes for the hotter 
stars. Also, it is clear from Figure [12] that the metal-rich stars lie 
somewhat low, possibly suggesting that the zero-point of the nii 
index as determined from Vega may be somewhat in error. 

When studying the calculated indices of Lester, Gray & 
Kurucz (1986) in Figure [TT] we find that they reproduce the 
observed sensitivity of the mi index more successfully than 
ours. However, these authors have transformed their calculated 
colours to match a set of standard stars and have thus scaled the 
amplitude of the mi index by a correction factor to fit the obser- 
vations. This is the probable reason for the closer agreement of 
those calculations with observations. We note that our line list 
is more complete, and that our treatment of the hydrogen line 
broadening (affecting H<5 and thus the v band) is more accurate. 



8 Onehag et al.: Calibration of Stromgren Piiotomtery 




I ■ ■ I ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ " ■ ■ I ■ ■ ■ ■ I I ■ ■ ■ I I ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I 




4500 5000 5500 6000 6500 7000 4500 5000 5500 6000 6500 7000 



eff ' eff 

Fig. 6. Theoretical Hp values versus effective temperature for models with different [Me/H] and \ogg plotted together with values 
for the standard stars (with [Fe/H] within +0.1 dex of the given [Me/H]). To guide the eye the theoretical log^ curves for 2.0 and 
4.5 are plotted as a dashed and a dotted line, respectively. 



Table 2. The metallicity sensitivity, 6(mi)/6[MefH]\i,^y), for 
models with log g = 4.5 & 4.0. 



(b-y)/lMem 


0.5-(-0.5) 


-l-(-2) 


-2-(-3) 


log = 4.5 








0.25 


0.124 


0.031 


0.008 


0.30 


0.145 


0.042 


0.013 


0.35 


0.149 


0.056 


0.028 


0.40 


0.139 


0.098 


0.026 


0.45 


0.136 


0.135 


0.032 


0.50 


0.176 


0.146 


0.048 


0.55 


0.233 


0.127 


0.070 


log ^ = 4.0 








0.25 


0.125 


0.030 


0.010 


0.30 


0.143 


0.041 


0.014 


0.35 


0.148 


0.059 


0.020 


0.40 


0.137 


0.081 


0.031 


0.45 


0.125 


0.097 


0.052 


0.50 


0.145 


0.100 


0.076 


0.55 


0.198 


0.093 


0.089 



What is then the reason for our discrepancy for early F stars? 
We have compared the MARCS model fluxes with observed so- 



lar and stellar fluxes from ground-based and space observations 
(Edvardsson 2008, Edvardsson et al. 2008) and traced probably 
significant discrepancies in the region 4000 A - 5000 A, with 
empirical fluxes of the Sun and solar-type stars being some- 
what smaller than model fluxes in the b band, while the blue- 
violet fluxes from HST/STIS of the more metal-poor stars are 
clearly in excess of the model fluxes in the violet-blue spectral 
region. These departures in both the b and v band may conspire 
to cause the discrepancy in calculated metallicity sensitivity. As 
discussed by Edvardsson et al. (2008), 3D model simulations 
suggest that these effects may be due to thermal inhomogeneities 
in the stellar atmospheres. Other systematic errors in the models, 
e.g. due to errors in opacities, line data and effects of departures 
from LTE are probably less significant. 

As a further test of our nii indices we will now be guided 
by the separate metallicity calibrations for F and G stars, re- 
spectively, by Schuster & Nissen (1989, eq. 2 & 3, hereafter 
S&N89). The derived F-star equation is based on the mi and 
(b - y) index and for the G-star equation the C[ index is also in- 
cluded. We derive a calibration expression based upon the form 
of these equations and using our theoretical colour grid. The 
results, when applying this to the standard stars and the C06 
sample, can be seen in Figures [T3] and [14] For F-stars we find 
standard deviations of derived [Me/H] values compared with lit- 
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Fig. 11. The mi index versus (b - y) for MARCS and Lester, Gray & Kurucz (1986, LGK) model atmospheres. Arrows show the 
effect of increasing the N abundances by a factor of 3. 



erature values for the standard-star sample as given in Table [T] 
We note, that for the standard sample, our theoretical calibration 
tends to suggest lower metallicities than the adopted values for 
a majority of the stars in the sample and increasing differences 
with increasing metallicity so that the differences amount to typ- 
ically 0.3 dex at [Me/H] = 0.0. The overall trend might indicate 
a zero-point problem; a shift of all stars by 0.1 30 dex gives a 
lower spread of the calculated metallicities compared to adapted 
values, (T+o = 0.217. This is however still not satisfactory in 
comparison with the empirical equations which generally show 
smaller spread and imply higher metallicities (see Table [T] and 
Figure [T3Tl. A linear regression for the F-star calibration values 
to the adopted values for the standard stars and the C06 sample 
is plotted in Figure [13] The standard deviation for this line fit is 
0" = 0. 172 dex. When applying our calibration to the C06 sample 
we see the same indications as for the standard sample, except 
for the lower metallicities, where the theoretical calibration sug- 
gest somewhat lower metallicities (-0.05 dex) than derived in 
C06. However, the 13 F-stars in this sample are too few to allow 
any definite conclusions. When making use of the 266 F-stars in 
the VF sample and remembering that this sample is highly biased 
towards metal -rich dwarf stars ([Fe/H]> -0.2 & \ogg > 4.0) we 
see the same trend as for the standard sample, i.e. our calibration 
suggests lower metallicities than those given by VF05 (see also 
TablelB- 



For the G-stars we find similar standard deviations for our 
standard sample with respect to literature values, see Table [T] 
There is an overall tendency that our theoretical calibration im- 
plies lower metallicities than listed in the literature for the metal- 
rich standard stars. The increasing deviation with increasing 
metallicity is of the same order as for the F-star calibration, i.e 
w 0.3 dex at [Me/H] =0.0. By adding 0.190 dex to all stars we 
would obtain a lower spread, cr+o 190 =0.181, which is of the 
same order as the deviation shown by empirical equations (see 
Table [1). A linear regression to the calculated metallicities for 
the standard sample and C06 is plotted as a solid line in the fig- 
ure (cr = 0.200). When using the 69 G-stars in the C06 sample, 
the result for the standard sample is confirmed, the theoretical 
calibration suggests lower metallicities for the higher metallic- 
ity range ([Fe/H] > -0.5) than listed in literature. Applying the 
calibrations to the 694 G-stars in the VF05 sample, we obtain 
similar results as for the C06 sample. 

For comparison the metallicities of our three comparison 
samples are calculated with 3 different empirical calibrations. 
The resulting standard deviations with respect to literature val- 
ues can be seen in Table [T] and linear regressions when applied 
to our standard sample can be seen in Figure[T3]&[l4]for F-stars 
and G-stars, respectively. 

We note in passing that Clem et al. (2004) also found dis- 
crepancies between their calculated mi indices and observations. 
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Fig. 12. The mi index versus (b - y) for MARCS models plotted together with standard stars of different gravities in different panels 
(logg for the stars within +0. 1 dex of the given log g). The metallicities of the stars are indicated by different symbols. 



For their coolest models they had to apply upward corrections to 
their calculated mi indices of 0.1-0.3 mag. Our nii indices de- 
part considerably less from observations but some upward cor- 
rections would be needed to fit the coolest stars in Figure fT2] 

The strong effects of microturbulence on the strengths of the 
dominating lines in the v band makes it important to investi- 
gate whether this could be the reason for the mis-match of the 
nil index in the theoretical calibration. In Table [3] we examine 
the effects of changing the microturbulence parameter The 
changes have been chosen to be larger for the lower-gravity mod- 
els to take the larger parameter values usually obtained by spec- 
troscopy for giants into account. We see that the effect of in- 
creasing the microturbulence is generally larger for lower log g 
and lower temperatures, as expected. When comparing to Figure 
[T2lwe see that this effect would indeed improve the situation by 
shifting the theoretical curves and steepening the curves in the 
low Teft end. Yet, the effects are considerably smaller than those 
needed to eliminate the mis-match. 

Another circumstance which might have some significance 
as an explanation for the problems with the mi sensitivity is the 
effect of lines from the (0,1) band of the CN violet system in the 
V band. In particular for the cooler giant stars, which may be af- 
fected by the first dredge-up of CNO-processed material, the CN 
lines may become stronger due to this; even if the carbon abun- 
dances are reduced by CNO processing, the enhanced N abun- 



Table 3. The change in the uvby indices when the microturbu- 
lence parameter is increased 



log ^ = 4.0 




= 1.0 




= 1.5) 




[Me/H] 


Aci 


Ami 


A{b-y) 


4500 


0.00 


0.0 


0.04 


0.011 


5500 


0.00 


0.003 


0.008 


0.005 


7000 


0.00 


0.004 


0.005 


0.001 


4500 


-1.00 


0.006 


0.005 


0.004 


5500 


-1.00 


0.006 


0.005 


0.002 


7000 


-1.00 


0.004 


0.002 


0.00 


log = 3.0 




A(e = 1.0 




= 2.0) 




[Me/H] 


Aci 


Ami 


A{b-y) 


4500 


0.00 


-0.016 


0.023 


0.032 


5500 


0.00 


0.004 


0.025 


0.014 


7000 


0.00 


0.001 


0.014 


0.002 


4500 


-1.00 


0.015 


0.018 


0.012 


5500 


-1.00 


0.017 


0.013 


0.003 


7000 


-1.00 


0.006 


0.005 


0.0 



dance (to which C is converted) makes the CN lines stronger. 
Also the '^CN lines should be significantly enhanced due to the 
production of '^C and N by CNO processing. We have explored 
these effects by systematically increasing the N abundances by a 
factor of 3, keeping the C abundance constant. This should lead 
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Fig. 7. Comparison of calibrations of (b - y) for dwarf stars with 
individual values for the standard stars and C06 plotted. The 
soUd line represents a fit of the effective temperatures adopted 
for the stars relative to the corresponding values obtained from 
the theoretical calibration. Below that (shaded area and right 
y-axis) corresponding linear regressions of the adopted effec- 
tive temperatures relative to empirical calibrations (Alonso et al. 
1996, Holmberg et al. 2007 and Ramirez & Melendez 2005) ai-e 
shown. The dotted line is a one-to-one line. 



Fig. 9. Comparison of calibrations of reff(H/?) for dwarf stars 
with individual values for the standard stars plotted. The solid 
line represents a fit of the effective temperatures adopted for the 
standard stars relative to the corresponding values obtained from 
the theoretical calibration. Below that (shaded area and right y- 
axis) a corresponding linear regression of the adopted effective 
temperautres relative to an empirical calibration (Alonso et al. 
1996) is shown. Note that the star located in the shaded area be- 
longs to the theoretical calibration and should therefore be read 
off at the left y-axis. The dotted line is a one-to-one line. 
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Fig. 8. Comparison of calibrations of (b - y) for giant stars with 
individual values for the standard stars plotted. The solid line 
represents a fit of the effective temperatures adopted for the stan- 
dard stars relative to the corresponding values obtained from the 
theoretical calibration. Below that (shaded area and right y-axis) 
corresponding linear regressions of the adopted effective tem- 
peratures relative to empirical calibrations (Alonso et al. 1999 
and Ramirez & Melendez 2005) are shown. The dotted line is a 
one-to-one line. 



to an overestimate of the effect, except for possibly stars high-up 
on the giant branch. As indicated in Figure [TT] this only leads to 
some effects for the more metal-rich giant stars and is not the 
explanation for the mis-match discussed here. 
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Fig. 10. The effective temperatures calculated by the equations 
based on H/3 and (b - y). Solid line theoretical fit, dashed line 
empirical fit (Alonso et al. 1996). The dotted line is the one-to- 
one relation. 



7. Surface-gravity calibration 

Theci = (M-v)-(v-y) index is designed to measure the Balmer 
discontinuity which is a temperature indicator for B- and A- 
type stars and a surface-gravity indicator for the late-type stars. 
Figure [15] shows its behaviour with changing parameters of the 
models. Obviously, it works nicely as a gravity criterion, with 
some dependence on metallicity for the most metal-rich stars, 
which has been a disputed issue in earlier calibration work. We 
note that for dwarfs cooler than the Sun it seems not very useful 
as a gravity measure, while for bright giants, and not the least 
metal-poor ones, it should work down to effective temperatures 
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Fig. 13. The [Me/H] calibration for F stars from the standard 
sample and Casagrande et al. (2006, C06). The solid line repre- 
sents a linear regression of the effective temperatures adopted for 
the stars relative to the corresponding values obtained from the 
theoretical calibration. Below that (shaded area and right y-axis) 
corresponding linear regressions of the adopted effective tem- 
perautres relative to empirical calibrations (Schuster & Nissen 
1989, Holmberg et al. 2007 and Ramirez & Melendez 2005) are 
shown. The dotted line is a one-to-one line. 
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Fig. 14. The [Me/H] calibration for G stars. For plot description 
see Figure [T3] Note that empirical calibrations of HNA07 and 
RM05 overlap and could therefore be difficult to distinguish. 



around 5000 K. In Figure [15] we have also plotted the indices 
calculated by Lester, Gray & Kurucz (1986). In view of the dif- 
ferences in line data and hydrogen-line theory we find the agree- 
ment satisfactory. 

Analogously with our treatment of the metallicity depen- 
dence of the nil index, we have measured the quantity 6ci /6 log g 
at constant {b - y) as a measure of the gravity sensitivity of ci . 
Again we have to turn back to earlier calibrations to find corre- 
sponding direct empirical measures. Thus, Schuster & Nissen 
(1989) have elaborated the methodology of Crawford (1975, 
1979) and write 



where 0.22 < (b - y)o < 0.47, (b - y)o being the dereddened 
(b - y), and 6co is the difference of a dereddened ci index and 
a standard sequence with Mv,zams at a given li/3. From this one 
may estimate the empirical sensitivity 6ci/6logg to be approxi- 
mately proportional to 2.5/ f. We thus obtain the following em- 
pirical values of{6ci/5logg; (b-y)) for [Fe/H] = 0: (0.21;0.30), 
(0.16;0.40) and (0.15;0.47). From Figure [15] and Table [4] we 
measure the corresponding theoretical 6ci/6logg values to be 
0. 19, 0.08 and 0.02, respectively. Thus, we find that the empirical 
gravity sensitivity of the ci index is well reproduced by the hotter 
models while it becomes underestimated for the cooler ones. We 
have also compared the calculated ci indices with the observed 
ones for our samples of standard stars and find a good agreement 
for the hotter stars (cf. Fig. [T6] l while for the cooler stars there 
is a severe mismatch. This is in itself not very remarkable - the 
total line blocking in the ultraviolet spectra of the cooler stars is 
considerably greater than 70%, and it is to be expected that this 
will not be very accurately described by the model spectra. The 
fluxes in the u and v bands are highly sensitive to other param- 
eters, such as microturbulence, and in particular CN abundance 
for the cooler giants, and these are known to vary systematically 
with gravity. In fact, as is indicated in Figure [15] the CN line 
strengthening which is expected for the red giants improves the 
fit to observed ci indices for the red giant models considerably. 

Clem et al. (2004) applied semiempirical corrections also to 
their calculated (u - v) and (v - b) colours in order to fit obser- 
vations. The resulting effect on their ci indices is typically less 
than 0. 1 mag while we would need a downward correction of our 
ci indices of about 0.2 mag for the coolest models to fit (cf. Fig. 
UM the metal poor stars and even more for the metal-rich ones. 



Table 4. The sensitivity to surface gravity, (6ci/dlogg)^i,-y), 
models with [Me/H] = 0.5 and 0.0 



for 



(b-y)/logg 


2.0-3.0 


3.0-4.0 


4.0-4.5 


[Me/H] =0.5 








0.25 


0.231 


0.216 


0.207 


0.30 


0.213 


0.187 


0.165 


0.35 


0.183 


0.148 


0.122 


0.40 


0.148 


0.112 


0.091 


0.45 


0.114 


0.086 


0.074 


0.50 


0.083 


0.071 


0.060 


0.55 


0.056 


0.062 


0.049 


[Me/H] =0.0 








0.25 


0.247 


0.230 


0.223 


0.30 


0.236 


0.205 


0.185 


0.35 


0.215 


0.169 


0.131 


0.40 


0.182 


0.126 


0.082 


0.45 


0.144 


0.081 


0.035 


0.50 


0.106 


0.040 


-0.006 


0.55 


0.065 


0.002 


-0.026 



Mv = Mv,ZAMS - f ■ Scq 

/ s 9.0 + 38.5 ■((/7-3;)o- 0.22), 



(2) 



In a recent article, Twarog et al. (2007) discuss the metal- 
licity dependence of the ci index for disk stars ([Fe/H] >-1.00) 
and conclude that the metallicity sensitivity of the index gener- 
ally has been underestimated. An interesting test would there- 
fore be to follow their recipe, i.e. divide our model grid into 
three groups (hot, warm and cold), and check whether or not our 
synthetic colours show the same behaviour For the hot mod- 
els ((b - y) < 0.43) we find the same strong metallicity de- 
pendence of the nil index and the same lack of sensitivity of 
the ci index as is found by Twarog et al. For warm models 
(0.43 < (b-y) < 0.50) we do find, like Twarog et al., a strong 
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ci sensitivity to metallicity, but likewise a relatively strong de- 
pendence of the nil index in contrast to the result of Twarog 
et al. For cool stars ((b - y) > 0.50), we separate unevolved 
dwarfs from subgiants and giants by making use of the defined 
LC (= ci - l.Omi + 3.0(b-y) - 0.15) calibration in Twarog et al. 
The models follow the same trends and LC separates the more 
metal-poor models ([Fe/H] =-0.5 & -1.0) but does not succed 
in separating the metal-rich ones ([Fe/H] > 0.0), as is also found 
by Twarog et al. It should however be noted that when apply- 
ing the metallicity calibrations for hot and warm stars as defined 
by Twarog et al. on our standard sample, we do not reach any 
higher accuracy in reproducing the literature values of [Fe/H] 
than when using e.g. S&N89. 

7.1. Applications 

In the area of application of the uvby photometry, some particu- 
lar questions have been of special interest to us: to which extent 
can the photometry be used for determining gravities and tem- 
peratures for metal-poor stars, and how sensitive is the synthetic 
ci index to certain elements such as nitrogen. 

Following Clem et al. (2004), we shall test our model colours 
versus observed colours for globular clusters. Unlike the ap- 
proach of Clem et al., our main emphasis lies on exploring the 



possibilities and shortcomings of our theoretical model colours 
sooner than correcting them semiempirically to establish new 
calibrations. 

7.1.1. M92 

As a check on the capabilities of the ci index for determining 
gravities, we have tested the model indices relative to the ob- 
served uvby photometry for the extreme Pop II globular clus- 
ter M92. Adopting a metallicity of [Fe/H]= -2.22, a distance of 
8.3 kpc, and a reddening of E{B - V) = 0.023 (Grundahl et al. 
2000), we have converted the observed y - (b - y) diagram of 
Grundahl et al. (private communication) to the fundamental pa- 
rameters along the evolutionary sequence past the turnoff point, 
using My, (b - y), Tss (determined by MARCS models) and as- 
suming a mass of O.8M0 for the giants; c.f. the Tf^g/ log g values 
given in Table |5] The bolometric corrections used were taken 
from VandenBerg & Clem (2003). Next, we have calculated val- 
ues for the uvby colours using MARCS models for these param- 
eter values, to compare with the directly observed dereddened 
indices for M92 stars, according to Grundahl et al. As is seen in 
Table|5] the observed colours are well reproduced by the models. 
We have also calculated the derivatives 6ci/6logg along the gi- 
ant branch and find typical values ranging from -0.2 to -0. 1 , thus 
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Fig. 16. The ci index versus (b-y) for MARCS models with different [Me/H] and log g plotted together with values for standard stars 
and Casagrande et al. (2006, C06) stars (marked as asterisks). The [Fe/H] intervals for the standard and C06 stars are [Fe/H] ±0. 1 dex 
and [Fe/H] +0.2 dex, respectively. 



agreeing within 20 % between the observations and the calcula- 
tions. This suggests that we can use the uvby system to estimate 
gravities for metal-poor giant stars. 



Table 5. Observed and calculated colours for logg values along 
the sub-giant - giant branch for M92 and MARCS, respectively. 
The tabulated colours are mean colour values of stars close to 
the given fundamental parameters. 









M92 




MARCS 








(b-y)o 


Co 


mo 




mi 


1.58 


4682 


0.584 


0.516 


0.109 


0.519 


0.117 


2.05 


4915 


0.522 


0.414 


0.088 


0.440 


0.088 


2.51 


5125 


0.473 


0.374 


0.075 


0.384 


0.067 


2.99 


5306 


0.437 


0.307 


0.069 


0.328 


0.052 


3.50 


5525 


0.400 


0.281 


0.057 


0.295 


0.047 


3.72 


5836 


0.355 


0.332 


0.047 


0.335 


0.045 


4.00 


6302 


0.299 


0.383 


0.064 


0.411 


0.055 



7.1.2. NGC6397 

As a second test, effective temperatures for stars in the metal- 
poor globular cluster NGC 6397 were derived. This is an ap- 
plication of current interest as Korn et al. (2007) have re- 
cently claimed that they have identified abundance trends be- 
tween groups of cluster stars which they suggest are caused by 
atomic diffusion. Tracing abundance differences between differ- 
ent groups, at the turnoff, on the subgiant branch and on the 
lower and higher red-giant branch (TOP, SGB, IRGB, hRGB), 
requires well determined stellar parameters, and most impor- 
tantly well measured effective temperature differences between 
the groups. Korn et al. derive effective temperatures with a 
photometric and spectroscopic approach and they find consis- 
tent temperatures. One potential problem, however, is the fact 
that two different photometric calibrations needed to be applied 
(Alonso et al. 1996, 1999) which may introduce systematic er- 
rors. For our theoretical calibration we derive effective temper- 
atures for each of the 18 stars in the four different groups, by 
using logg values and dereddened colour indices {b - y) and 
(v - y) listed in Korn et al. and adopting their [Fe/H] of -2.0. 
Mean values of the effective temperatures for the four groups are 
shown in Table |6] The results indicate that our effective temper- 
atures are 50K to 1 00 K hotter than obtained from the empirical 



dnehag et al.: Calibration of Stromgren Photomtery 



15 



Table 6. Derived effective temperatures from empirical calibra- 
tions of photometry {{b - y) and ( v - y)), and spectroscopy (Korn 
et al. 2007), as well as our theoretical calibration, for groups of 
stars in the metal -poor cluster NGC 6397. The temperatures are 
mean values for a number of stars/models in each group. 





Phot. 




Spec. 


MARCS 




Teff [K] 


(h-y) 


(V-.V) 




(b-y) 


(v-y) 


hRGB 


5121 


5132 


5130 


5225 


5226 


IRGB 


5455 


5408 


5456 


5541 


5535 


SGB 


5797 


5824 


5805 


5839 


5853 


TOP 


6229 


6214 


6254 


6288 


6281 


A(TOP-RBG) 


1108 


1081 


1124 


1063 


1055 


A(TOP-SGB) 


433 


390 


449 


449 


428 



Table 7. Calculated 6cy/6[NfFe] for giant models as a compar- 
ison to the observed Cy sensitivity of NGC 6752 (Yong et al. 
(2008)) 





log,? 


0.0—0.75 


0.75—1.5 


4749 


1.95 


0.03 


0.00 


4829 


2.10 


0.04 


0.02 


4841 


2.15 


0.04 


0.02 


4904 


2.33 


0.04 


0.04 


4950 


2.42 


0.04 


0.04 



calibrations. More importantly, the temperature differences be- 
tween TOP, SGB and hRGB stars are close to (even somewhat 
smaller than) the ones obtained from the empirical calibrations 
by Korn et al. This strengthens the claim of Korn et al. that abun- 
dance differences in Li, Mg and Fe, reflecting atomic diffusion, 
are present at a significant level. 

7.1.3. NGC 6752 

In another test, the nitrogen sensitivity of of the ci index was 
examined. In a recent paper, Yong et al. (2008), measure the 
nitrogen content of giants in the globular cluster NGC 6752 
([Fe/H] ~ -1.6) based on high-resolution observations of NH. 
The authors find a strong correlation between nitrogen content 
and the ci index and establish the c,, index (ci - (b - y)) that 
shows a close to linear correlation with respect to nitrogen con- 
tent. To examine the synthetic sensitivity of this c,, index, a set of 
models and high resolution spectra were calculated (Teff = 4749, 
4829, 4841, 4904, 4950 K, log^ = 1.95, 2.10, 2.15, 2.33, 2.42, 
[Me/H] - -1.50) based on a selection of stars discussed by Yong 
et al. The [N/Fe] content was set to vary as 0.00, 0.75 and 1.50 
for each of the models in order to cover a great part of the [N/Fe] 
span in NGC 6752 found by Yong et al. As a result we find that 
our synthetic Cy index is sensitive to the nitrogen content of the 
model, even though the sensitivity 6cy/6[N/Fe] of the models is 
somewhat smaller, ~ 0.03 - 0.04 (see Table|2ll, compared to the 
sensitivity Yong et al. find for their giant stars, ~0.06, as is esti- 
mated from the best-fit straight line in Yong et al. Figure 7. 

8. Conclusions and recommendations 

We have explored the possibilities and shortcomings of syn- 
thetic uvby-H/S photometry based on new MARCS model atmo- 
spheres. In general a good agreement with empirical calibrations 
of this photometric system is found. However, a number of sys- 
tematic deviations between theory and observations also become 
apparent. The temperature sensitivity of the (b - y) colour (i.e. 



6(b -y)l dTeff) seems marginally larger for the calculated colours 
than is found when using infrared-flux-method determinations 
of temperatures. A similar, and even somewhat larger, difference 
occurs for temperatures based on the H/3 index, when compared 
with the empirical scale. The (b - y) calibration for Pop II stars 
supports Korn et al. (2007), who claim the signatures of atomic 
diffusion in the metal-poor globular cluster NGC 6397. 

For the metallicity index mi the theoretical sensitivity 
6mi/6[Me/H] is somewhat larger than the empirical one that 
is based on spectroscopic [Fe/H] -determinations. For the grav- 
ity sensitivity of the ci B aimer discontinuity measure, we find 
a reasonably good agreement with observations for stars hot- 
ter than the Sun, where the Balmer discontinuity is significant. 
Considerable problems remain for the cooler stars, although the 
model calibration works well for Pop II giants. 

One might ask whether these problems may be solved when 
even more detailed atomic and molecular line data become avail- 
able to feed into model atmospheres. This is possible but it is 
also possible that thermal inhomogenities as generated by con- 
vection and possibly also non-LTE effects, contribute to the dif- 
ference between synthetic and observed values of (b-y) and mi . 
Concerning the problems in reproducing the observed ci indices 
of cooler stars, these may well be due to failures in line-data, 
although thermal inhomogeneities in the atmospheres may also 
here turn out to be the important effect. One should also real- 
ize that there are still residual problems associated with differ- 
ences and uncertainties in the observed and adopted Stromgren 
passbands (Manfroid & Sterken, 1987). The current study illus- 
trates the problems in synthetic photometry in the visual wave- 
length regions. Continued efforts along such lines must be com- 
plemented with more detailed studies of the shortcoming of clas- 
sical model atmospheres and the replacement of those, also in 
large-scale calibration efforts, by physically more realistic mod- 
els. 
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Appendix A: Calibrations 

In the following Tcs and [Pe/H] calibrations we have systematically used 
the form, and in most cases the index Hmits for the validity of the calibra- 
tions, according to Alonso et al. (1996,1999) and Schuster & Nissen (1989), 
respectively. The coefficients have then been calculated by a least square 
method to optimize the fit to the model indices. 



A.1. Theoretical Teff calibration for dwarf stars, (b - y) and 

Oeff = 5040/Teff 

0eff(/'->')dwarf = 0.466 + 1 .025(* - y) - 0. 123(Z) - j;)^ 
-1-0.212(6 - y)ci - 0.051(6 - >')[Me/H] 
-(■0.003 [Me/H] - 0.005 [Me/H]^ 

The models, upon which the calibration is based, are selected to follow 
the restrictions given for the empirical equation (A96, Eq. 9. Restrictions 
adopted from Rig. 11a): -3.00 < [Me/H] < 0.50; 0.70 < {b - y) < 0.20; 
0.10 < ci < 0.55; 4.00 < logg (set as limit for dwarf stars). 



©etf (HyS)dwaif = 28.60 - 19.79H/3 + 3.504H/32 + 0.422HyS[Me/H] 
-1.068[Me/H] + 0.002[Me/H]^ 

In analogy with the theoretical calibration based on (b - y), the limits 
determining which models to use, were set to be equal to the ones for 
empirical calibration (A96, Eq. 10 and its following applicable ranges): 
2.44 <Hy3 < 2.74 for 4-0.5 > [Me/H] > -0.5; 2.50 <H;3 < 2.70 for 
-0.5 >[Me/H]> -1.5; 2.50 <H/3 < 2.63 for -1.5 > [Me/H] > -2.5; 
2.51 <H/3 < 2.62 for -2.5 > [Me/H] > -3.5; 4.0 < logg (set as a limit 
for dwarf stars). 



A.2. Theoretical Teif calibration for giants, (b - y) 

©cff (^ - 3')giant = ao + ai{h-y) + a2{b-y)- + 

03(6 - )')[Me/H] + a4[Me/H] + a5[Me/H]^ 



Table A.l. The derived coefficients for the theoretical calibration 
in the intervalls I: 0.15 < (b - y) < 0.424, II: 0.424 < (b - y) < 
0.712 & -5.0 < [Me/H] < -0.5 and IH: 0.428 <ib-y)< 0.794 
& -0.5 < [Me/H] < 0.5 





ao 


ai 


32 


as 


a4 


as 


I: 


0.6732 


0.0859 


1.1455 


-1.080e-2 


-0.132e-2 


-0.082e-2 


11: 


0.1983 


2.0931 


-0.9978 


4.709e-2 


-2.66e-2 


-O.Ole-2 


III: 


0.4522 


1.1745 


-0.3093 


-0.1693 


2.165e-2 


-1.679e-2 



Here we deviate from the same restriction limits for the calibrations as pre- 
sented in Alonso et al (1999, Table 2 and 3). New theoretical limits were set 
due to model restrictions: (6 - y) < 0.424 and (h — y) > 0.424 (cf. Alonso 
et al 1999, Eq. 14: (6 - y) < 0.55 and Eq. 15: (6 - y) > 0.50). The logg 
limit was set to 1.5 < logg < 3.5. 



A.3. Theoretical metallicity calibration for G-stars 

The models, upon which the metalHcity calibrations for G and R stars 
are based, are selected to follow the restrictions given for the empirical 
equations of S&N89, Eq. 3. and 2., respectively. The G-stars limits cover 
the following intervals: -2.6 < [Me/H] < 0.4; 0.37 < (6 - y) < 0.59; 
0.03 < ci < 0.57; 3.00 < log g (set as limit for log g). 



[Me/H]G = -2.796 -t- 39.21mi - 88.97m^ - 73.43mi(6 -y) 

-i-181.4m^(6-y) -i- [27.03mi - 1.220ci -41.42mJ]ci 
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A.4. Theoretical metallicity calibration for F-stars 

F-star limits are given by: -3.5 < [Me/H] < 0.2; 0.22 < (b - y) < 0.38; 
0.17 < ci < 0.58; 3.00 < logg (set as limit for logg). 

[Me/H]F = 1.850-34.21mi + 105.34mi(/)->') 
+ n9.8mj{b -y)- 242.4mi(/> - yf 
+[2.757 - 20.38mi + Q.llllib - y)} log(mi - C3) 

C3 = 0.4462 -2.233(fe-y) + 2.885(fe-y)^ 



